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Scanning SQUID susceptometry images the local magnetization and susceptibility of a sample. 
By accurately modeling the SQUID signal we can determine physical properties such as the pen- 
etration depth and permeability of superconducting samples. We calculate the scanning SQUID 
susceptometry signal for a superconducting slab of arbitrary thickness with isotropic London pen- 
etration depth A, on a non-superconducting substrate, where both slab and substrate can have a 
paramagnetic response that is linear in the applied field. We derive analytical approximations to 
our general expression in a number of limits. Using our results, we fit experimental susceptibil- 
ity data as a function of the sample-sensor spacing for three samples: 1) 5-doped SrTiOa, which 
has a predominantly diamagnetic response, 2) a thin film of LaNiOs, which has a predominantly 
paramagnetic response, and 3) the two-dimensional electron layer (2-DEL) at a SrTiOa/LaAlOa 
interface, which exhibits both types of response. These formulas will allow the determination of 
the concentrations of paramagnetic spins and superconducting carriers from fits to scanning SQUID 
susceptibility measurements. 

PACS numbers: 74.72. Cj 85.25.Dq 74.25.Ha 



I. INTRODUCTION 

Scanning SQUID microscopy^^ allows the simultane- 
ous imaging of the local magnetization and the magnetic 
response (susceptibility)^ of the surface of a sample on 
a micron length scale. The sign and magnitude of the 
susceptibility signal yields information about electrons 
in the material. For a superconductor, the diamagnetic 
susceptibility is a measure of the local London pene- 
tration depth.'^K^ In most superconductors, the diamag- 
netic susceptibility is much stronger than other possible 
sources of magnetic response, such as nuclear suscepti- 
bility or the paramagnetism of impurities, other regions 
of the sample, or non-superconducting carriers. How- 
ever, in superconductors with unusually strong compet- 
ing paramagnetic susceptibility and/or a low superfluid 
density, it may be necessary to consider both types of 
contributions. For example, a paramagnetic response 
has been observed in scanning susceptometry measure- 
ments of non-superconducting sampled and supercon- 
ducting samples above their critical temperatures.^^^ 

The temperature dependence of the London penetra- 
tion depth, which is related to the susceptibility, has 
played an important role in determining the symmetry of 
the superconducting order parameter in unconventional 
superconductors.^^ However, for superconductors with 
low superfluid densities, the diamagnetic contributions 
from Cooper pairs and the paramagnetic contributions 
from spin or other sources can have similar magnitudes 
but different temperature dependences, making it diffi- 
cult to determine the temperature dependence of the su- 
perfluid density. It is therefore important to be able 
to separate the paramagnetic from the superconducting 
components in scanning SQUID susceptometry measure- 



ments. 

Kogan^ presented a model for the diamagnetic re- 
sponse of a superconductor to arbitrary local field 
sources. One source he considered was a circular ring 
of current appropriate for scanning SQUID susceptome- 
try. Here we extend his model to include both diamag- 
netic and paramagnetic effects, for a paramagnetic super- 
conductor of arbitrary thickness on a paramagnetic sub- 
strate. Our final expression reduces to that of KogarP^ 
for a superconductor with the permeability of vacuum in 
the bulk and thin film limits, and to that of Bluhm et 
alW for a thin film paramagnetic response. We present 
in Table [l| analytical approximations to our full expres- 
sion for a) a bulk non-superconducting paramagnet, b) 
a thin non-superconducting paramagnet, c) a bulk su- 
perconductor without paramagnetism with penetration 
depth short relative to the other lengths in the problem, 
d) a bulk superconductor without paramagnetism with 
penetration depth long relative to other lengths in the 
problem, and e) a thin superconductor without param- 
agnetism. These analytical approximations, along with 
the full expression Eq. [7| represent the main results of 
this paper. 

Although in this paper we concentrate on the scan- 
ning SQUID susceptometer geometry, with the field coil 
co-planar and co-axial with the pickup loop , the same 
basic formalism could be applied to penetration depth 
measurements using other two-coil mutual inductance 
geometries ii^tti^J fQj. example, with the two coils at dif- 
ferent heights, or on opposite sides of a thin film sample. 

As examples of applications of these expressions we 
fit scanning susceptometry data on a ^-doped sample of 
SrTiOs, a thin film of LaNiOa, and the two-dimensional 
electron layer (2-DEL) at the interface between SrTiOa 
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FIG. 1: Model geometry. a) The layout of the field 
coil/pickup loop region of the SQUID susceptometers used 
in this paper, b) Approximations to this layout considered 



in Appendices 1 a and lb c) For this paper we consider 
a slab of thickness t with magnetic permeability /J2 and su- 
perconducting penetration depth A on a non-superconducting 
substrate with magnetic permeability /is- We define ^ = at 
the sample surface, with a pickup loop (radius b) and concen- 
tric field coil (radius a) located at a height z = zq in a plane 
parallel to the sample surface. 



and LaAlOsPFor the case of (5-doped SrTiOa, which has 
a predominantly diamagnetic response, there is a strong 
correlation between the best penetration depth and sen- 
sor height parameters in fitting susceptibility approach 
curves: an uncertainty in the sensor height of 1.5/im re- 
sults in an uncertainty in the Pearl penetration depth 
A = 2A^/t of 60%. A susceptibility approach curve for 
the LAO/STO interface at a position which is predom- 
inantly paramagnetic is well fit by our expression for a 
thin film paramagnet with reasonable values for the fit- 
ting parameters. However, at a position which is weakly 
diamagnetic the best fits for the height and permeability 
parameters take on unphysical (negative) values, even for 
a model which includes both superconducting and para- 
magnetic contributions. We speculate that this last may 
be due to sample inhomogeneity and/or an interaction 
between the SQUID and the sample superfluid. 



redirection from the leads.^In this paper we assume the 
geometry of Fig. [1J3: the susceptometer is represented 
by two co-planar concentric circular loops. The field coil 
has radius a, and the pickup loop has radius b . Both 
are infinitely thin wires. We evaluate in Appendices 1 1 a| 
and | 1 b| the systematic errors in the SQUID susceptibility 
associated with our approximations to the field coil and 
pickup loop shapes. The loops are oriented parallel to, 
and a height zq above, a slab of material with thickness 
permeability /i2 and isotropic London penetration depth 
A on a semi-infinite non-superconducting substrate with 
permeability jus- 

One can divide space into 3 regions of interest: 1: 
zo > z > 0, 2: > z > -t, and 3: -t > z. The 
half-space —t>z has permeability jus and is non super- 
conducting. The space t > z > has /ii = jno and is also 
not superconducting. 

Following Kogan,^^ V x = in region 1 since there 
are no currents. This identity is assured by writing the 
field Hi as the gradient of a scalar potential cpi: Hi = 
S/ifi. Then V^i = since V • ^1 = V • (/io^i) = 0. In 
region 3 there is again no transport current, so V x Hs = 
0. Since we consider only the case where ja is spatially 
homogeneous, V • Hs = 0, and therefore \/'^(f3 = 0. 

Inside the superconductor the total current, propor- 
tional to the curl of the magnetic flux density consists 
of the supercurrent plus the current due to the inhomo- 
geneity of the magnetization M: 

where the London penetration depth A is assumed 
isotropic and homogeneous, is the quantum mechanical 
phase, A is the magnetic vector potential, and ^0 = h/2e 
is the superconducting flux quantum. Taking the curl, 
the magnetic flux density B satisfies^^ 



V X V X ^ + B/X^ /io V X V X M. 



(2) 



We here neglect a term which represents a sum of delta 
function vortices. If we consider the case where there 
is a small susceptibility such that M = and write 
= M2//^o = 1 + X2, with X2 << 1, we then recover the 
familiar London's equation 



II. MODEL 



A. Full expression 



We consider the geometry of Figure [T] The SQUID 
susceptometers used in this paper have the layout shown 
in Fig. [1^.^^ Traditionally such a layout has been approx- 
imated by that of Fig. [T]3: the field coil is represented by 
a circular wire, while the pickup loop is represented by a 
circular wire plus an additional pickup area due to flux 



X^V^B -5 = 0, 

with a modified penetration depth A = X/^/p^ 
all the experiments reported here |/i2 — 1| << 



(3) 



1: the 



i) are less than 
3le for the para- 



fit values for X2^ (see, e.g. Fig. |13^ 
6xl0~^/im. Even if the layer responsi : 
magnetism is only lOnm thick, this would correspond to 
X2=p2-l< 0.06. 

The fields in the 3 spatial regions of interest can then 
be expanded in Fourier series as 



(pi{f,z) 



ik-' 
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^2(r,z)- ^ 



(4) 



where if>s is the source potential due to currents in the 
susceptometer field coil, and (^r,ii V^r,3, and H~ 
are response potentials and fields, k = |^|, and q = 



(k'^ + A~^)^/^. Applying the boundary conditions of con- 
tinuity of the normal component of B and the tangential 
component of H at the interfaces 2; = and z = — t, as 
well as the requirement that V • 5 = in region 2, leads 
to the solution 



ifr,l{k) 



-{q + k/j.2){-kfi2 + qfis) + - k/j.2){k/j.2 + q/J^s) 

-{q - kfl2){-kfl2 + qp^s) + e'^'i^{q + kfl2){kfl2 + ^^3) 



(Ps{k) 



(5) 



Simplified versions of this expression are given in Table 
[l|in five limiting cases: a) is a bulk, non-superconducting 
paramagnet, b) is a thin, non-superconducting paramag- 
net, c) is a bulk superconductor without paramagnetism, 
with penetration depth short relative to the sensor field 
coil radius and height, d) is a bulk superconductor with- 
out paramagnetism, with penetration depth long relative 
to the field coil radius and height, and e) is a thin super- 
conductor without paramagnetism. It is of interest to 
note that in three of the cases: bulk paramagnetic (a), 
thin paramagnetic (b), and bulk weak diamagnetic (d), 
the material property of interest, either the permeability 
/i2 or the penetration depth A, is separable from a ge- 
ometrical factor, independent of the form of the source 
potential. This means that the temperature dependence 
of the material property can be determined, aside from 
a multiplicative constant, without curve fitting, in these 
cases. For example, in the bulk weak diamagnetic case 
(Table llli), the SQUID susceptibility is proportional to 



A~^, with a constant of proportionality that depends only 
on geometry, independent of the form of the source term, 
which should be independent of temperature. This is 
also true in the thin diamagnetic case (Table |l^) for suf- 
ficiently large Pearl lengths A. A and the geometrical fac- 
tors are not separable in the limiting case of bulk strong 
diamagnetism (Table [i fc). I n this case, to a good approx- 
imation (see Appendix 1 1 c[ ) it is the sum A + zq which is 
determined by SQUID susceptibility measurement s.^^^ 

The source field for a circular field coil of radius a is 
given by^ 



(fs{k) 



irla 



-^Vi(te) 



(6) 



The 2;-component of the response field in region 1 is given 
by hr{k^z) = —k(pr^ie~^^ . Taking the limit b « a, the 
height dependence of the SQUID susceptibility (j){z) is 
given by 



-f 

Jo 



dxxe 



-2xz 



Ji{x) 



-{q - fi2x){lJ.3q - l^2x) + e^^*{q + iJ.2x){iJ.3q + IJ.2x) 



(7) 



where 



1 

^~dl' 



^ is the fiux through the pickup loop in response to 
the current /, ^0 = h/2e is the superconducting fiux 
quantum, the self inductance between the field coil and 
the pickup loop 



A/io/2<^>oa, 



A is the effective area of the pickup loop, z = z/a, t 
t/a^ and q = \/x^^^Y/)^ ^ with A = A/a. 



B. Analytical expressions in various limits 

(8) In general the integral in Eq. [7|must be done numer- 
ically, but analytical expressions can be derived in the 
limits given in Table |T] and plotted in Figure |2j 

It is to be noted that in the bulk limit {t » 1) the 
strong superconducting susceptibility (A << 1, /i2 = 1) 
and the paramagnetic susceptibility (A >> 1, fl2 > 1) 
have the same height dependence aside from a scaling 
factor — (/i2 — 1)/(M2 + !)• However, when the thickness 

(9) of the paramagnetic superconductor becomes compara- 
ble to the field coil radius, the height dependences of the 
paramagnetic and superconducting contributions become 
different, and it is possible in principle to determine the 
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Description Thickness Penetration depth Permeabihty (pr,i{k) /(ps{k) 



Bulk para 
Thin para 
Bulk strong dia 
Bulk weak dia 
Thin dia 



t»l,z 
t «l,z 

t » 1 
t» l,z 
i «l,z 



A >> l,z 
A << l,z 
A >> l,z 
X»t 



M2 > 1 

/i2 > 1, As = 1 

A2 = 1 

A2 = As = 1 
A2 = As = 1 



A2+1 
A2 + 1 
fct(/x2_i) 

1 

A2(g+fc)2 
1 

4A2fc2 

1 

1+Afc 



VA2 + 1/ (1+422)3/2 

Aj-l 3tz 
A2 (1+422)5/2 

1 

(l+ 4(2 + A)2) 3/2 

^(74^^-2^) 



-f (l 



Vl+42 



TABLE I: Response scalar potential (pr,i{k) divided by the source potential (ps{k) in momentum space (Eq. |5| and scanning 
SQUID susceptibility (/)(z) divided by the SQUID self-susceptibility (/)s in real space (Eq. [7|, in various limits. 



a) Bulk paramagnetic 

b) Thin paramagnetic 

c) Bulk strong diamagnetic 

d) Bulk weak diamagnetic 

e) Thin diamagnetic 




a) 



6-doped STO 



3 b) x10"^ 



FIG. 2: Theoretical height dependence of the scanning 
SQUID susceptibility, divided by the maximum of the ab- 
solute value of the susceptibility, of a paramagnetic supercon- 
ductor in various limits. The letters correspond to the entries 
in Table m 



relative contributions of each to the total susceptibility 
by fitting approach curves. It is also possible in princi- 
ple to determine the z-dependence of the response car- 
rier density (either paramagnetic or diamagnetic) from 
approach curves. However, in practice the differences 
between the spacing dependences of the various contri- 
butions are subtle, and it is difficult to separate out the 
paramagnetic from the diamagnetic components without 
extra information. Such information could be supplied, 
for example, by raising the temperature above the super- 
conducting transition temperature, leaving only the para- 
magnetic contribution, or studying the low temperature 
dependence of the susceptibility, where the temperature 
dependence of the superconducting component could sat- 
urate, while that of the paramagnetic component could 
become larger. Finally, one or both components could be 
spatially dependent (see e.g. Fig. |4|, which could help 
to separate them. The regions in parameter space of va- 
lidity and errors associated with using the approximate 
expressions in Table [l| are explored in Appendix 1 1 c[ 
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FIG. 3: a) Susceptometry image of a ^-doped SrTiOs sample, 
b) Susceptibility as a function of AVz, the change in z-piezo 
voltage from contact between the SQUID substrate and the 
sample surface, at the position of the gray circle in a). The 
dots are data, the line is a fit using the thin diamagnetic 
limit expression in Table IJe), with A = 954/im. c) scan- 
ning susceptibility image of a patterned non-superconducting 
thin film of LaNiOs. d) Susceptibility approach curve for the 
LaNiOs film at the position of the square symbol in (c). The 
dots are data, the line is a fit using the thin paramagnetic 
limit expression (Table [l];b)), with X2t = 1.3 X 10 fim. 



III. COMPARISON WITH EXPERIMENTS 

There have been a number of works in which SQUID 
susceptibility measurements have been used to infer 
the London penetration depth of superconductors .I^H^ 
We examine scanning SQUID data from several sam- 
ples. The low temperature measurements were per- 
formed in a home-built SQUID microscope in a dilution 
refrigerator.!^ The 5K measurements were performed in a 
home built variable sample temperature scanning SQUID 
microscope.^ The SQUID susceptometers used in both 
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microscopes were described in Ref. [T8l 

Figure |3] shows experimental data for samples with pre- 
dominantly diamagnetic response (Fig. IsLb) and para- 
magnetic response (Fig. |3]3,d). Fig. ^3] (a) and (b) 
show SQUID susceptometry of a Nb (5-doped sample of 
SrTi03(ST0).'^This sample was grown in an atmosphere 
of 10~^ Torr oxygen at 1200° C. Nb dopants were con- 
fined to a 5.9nm layer, with 100 nm cap and buffer lay- 
ers of STO grown above and below the doped region. 
The sample was annealed in situ under an oxygen par- 
tial pressure of lO-^Torr at 900° C for 30 minutes.^^ For 
the data sets of Fig. [s] (b) and (d), the susceptibility 
was recorded while the SQUID was driven towards the 
sample by increasing the z-piezo voltage Vz from a large 
negative value. In these plots /SVz = corresponds to 
contact between the SQUID substrate and the sample 
surface. At contact the SQUID pickup loop is a height 
zq above the sample surface because of the finite angle 
(typically a few degrees) between the sample surface and 
the SQUID substrate surface. The dots in Fig. [sj^b) 
show such a susceptometry approach curve at the posi- 
tion indicated by the gray circle in Fig. [sj^a). The line 
(difficult to distinguish from the data in this plot) is a fit 
of this data to the thin diamagnetic limit of Table |T] with 
5 fitting parameters: a vertical shift 3(j)^ a linear slope 
^linear = ^LZ^ the Pearl length A, zq, and the change in z 
with piezo voltage dz/dVz- We do not know the source of 
the linear background. In the present case it was small, 
a - -3.7 X 10"^ l/mA-/im. 

The two fixed parameters in this analysis were the ef- 
fective field coil radius a=8.4/im, and pickup loop radius 
6=2.7/im. The effective field coil radius was taken from 
the numerical calculations of Brandt and Clem,^^ using 
(see Fig. [l^) a field coil inside radius of 6.5/im, outside 
radius of 12/im, thickness 0.3/im, and penetration depth 
0.09/im. The effective pickup loop radius was chosen to 
result in the measured self inductance of (j)s = 800 1/A 
using Eq. |9] This results in an effective pickup loop area 
of 22/im^, larger than the 17/im^ obtained from the sum 
of the geometric mean of the pickup loop itself, with in- 
side radius ri^=0.88/im, and outside radius ront=2.4/im, 
added to the Ketchen's 1/3 rule area (w^/3)[T9] for the 
shield over the pickup loop leads, which has w;=4.5/im. 
Part of this discrepancy may be due to the fact that the 
pickup loop shield focusses flux from the fleld coil into 
the pickup loop area. 

The paramagnetic susceptibility of samples is typically 
much smaller than the diamagnetic susceptibility of su- 
perconductors. An example is shown in Fig. [sjc) and 
(d). Fig. [Sjc) shows a scanning susceptometry image of 
a patterned 20nm thick fllm of LaNiOa, imaged at 5K 
with a SQUID susceptometer with the same geometry 
as in Fig. [sja). For this sample thermally evaporated 
La and Ni were co-deposited on to a LaAlOa substrate 
kept at a temperature of 600C, in a background oxy- 
gen pressure of 7 x 10~^ Torr. An RF source operated at 
200 W provided atomic oxygen to the fllm during growth, 
and in situ structural characterization was obtained using 
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FIG. 4: Scanning SQUID susceptometry image of a patterned 
LAO/STO interface at 0.087K. The labels indicate where the 
data in Fig. [5] was taken. 
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FIG. 5: Scanning SQUID susceptibility as a function 
of sensor-sample spacing at 2 positions on the patterned 
LAO/STO sample of Fig. [3] PI, in a region showing pre- 
dominantly diamagnetic shielding; and P2, in a region show- 
ing predominantly paramagnetic response. The symbols are 
data and the lines are fits as described in the text. 



RHEED. Patterning of the fllm was achieved by creating 
a mask on the fllm surface using photolithography, and 
then etching the fllm in a HCl solution (4:1 H20:HC1) 
to remove uncovered areas. The dots in Fig. |3]i are 
the data, the line is a flt to the thin paramagnetic limit 
expression of Table |l]^b), with 4 fltting parameters: (50, 
X2t = 1.3 X 10~^/im, zo=1.9/im, and dz/dVz=2.7fim/V, 
where X2 = /mo — 1 and t is the thickness of the fllm. 

As an example of a sample that shows both paramag- 
netic and diamagnetic behavior we present data on the 
2-DEL at the interface between the perovskite insulators 
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SrTiOs (STO) and LaAlOs (LAO). The sample and mea- 
surement techniques for the data used in this study were 
described in Ref. 9, Briefly, the sample was prepared by 
growing 10 unit cells of LaAlOs on a commercial Ti02 
-terminated 001 STO substrate, with an aluminum ox- 
ide hard mask patterned on to the STO substrate prior 
to LAO growth. A crystalline LAO/STO interface only 
grew in the gaps of the patterned mask. The LaAlOs was 
deposited at 800° C with an oxygen partial pressure of 
10~^ mbar, after a pre-anneal at 950° C with an oxygen 
partial pressure of 5 x 10~^ mbar for 30 min. The sam- 
ple was cooled to 600° C and annealed in a high-pressure 
oxygen environment (0.4 bar) for one hour. Figure [s] 
displays SQUID susceptibility data as a function of spac- 
ing between the sensor and the sample surface for the 
LAO/STO sample imaged in Fig. [4] at the positions la- 
belled. Both positions are in a gap of the aluminum oxide 
mask, but P2, close to the edge of the two-dimensional 
electron layer, shows paramagnetic response, while PI 
shows diamagnetic behavior. These approach curve data 
were taken at T=0.02K with a field coil current of 1mA. 

The fact that P2 shows a maximum below AVz = 
implies that the paramagnetism results from a thin 
film, rather than from the substrate (compare the thin 
and bulk paramagnetic limit curves in Fig. |2|. Fit- 
ting this data to the pure thin paramagnetic expression 
of Table [ij with X2^, ^o, <^ and Scj) as variables, with 
a=8.4/im, 6=2.7/im and dz/dV = 2.9/im/V, results in 
X2^=4.9+0.8-0.7x 10-"^ /im and zq = 1.5 + 0.7 - 0.3/im. 
This fit is displayed as the solid line in Fig. |5j Using the 
same assumptions as for the LaNiOa case above (includ- 
ing a ±20% uncertainty in a), we find a spin density of 
1.25±0.5xl0^'^l/cm2. 

We attempted to fit curve PI in Figure |5] to the pure 
thin film diamagnetic expression of Table |l]^e), with A, zq, 
a, and S(j) as variables, and dz/dV = 2.9/im/V as a fixed 
parameter. The best fits were obtained for unphysical 
negative values for zq . If we constrain zq to vary between 
the values of 1 and 2.5/im, the best fit (dashed line in Fig. 
[5| occurs for zq = 2.5/im and A = 16.4/im. However, 
the fit quality was not good (the best fit S^-^ is about 
25 times worse for PI than for P2). Using the same 
assumptions as for ^-STO above, but with an effective 
mass m* = 1.46me^^, the allowed values for A (15 mm 
< A < 34mm) correspond to a Cooper pair density of 
lxlO^^cm-2 < Ns < 3.4 x lO^^cm-^ 

It seems reasonable to assume that the susceptibility 
at position PI in Fig. [4] has both superconducting and 
paramagnetic contributions, and therefore could be fit 
using the full expression Eq. [T] However, the fitting 
parameters X2^, and A are strongly correlated, result- 
ing in large uncertainties in their values. Further, fits to 
this data result in unphysical negative best fit values for 
X2t and Zq. We speculate that these unphysical values 
might result from the inhomogeneous superfiuid density 
in this sample or from interaction between the sensor 
SQUID and the superfiuid at these low densities. There- 
fore, as mentioned in the introduction, additional infor- 



mation, such as different temperature or spatial depen- 
dences, will be required to separate the superconducting 
from the paramagnetic components in scanning SQUID 
susceptibility measurements. 



IV. CONCLUSIONS 

We have presented a full expression and analytical ap- 
proximations in various limits for the susceptibility in 
a scanning SQUID geometry of a paramagnetic super- 
conductor of arbitrary thickness on a paramagnetic sub- 
strate. These expressions can be used to measure the 
spin concentration and the Cooper pair density in a para- 
magnetic superconductor. A comparison of analysis 
with bootstrap statistical analyses (see e.g. Appendix 
Tdl Fig. 



12) indicate that the accuracy of these mea- 



surements can be improved with a precise knowledge of 
the sensor height Zq and the piezo constant dz/dV in 
scanning SQUID susceptometry measurements. 
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Appendix: Sources of error 
1. Systematic 

a. Approximating field coil by circular wire 

The actual susceptometer layout used in the experi- 
ments described in this paper is shown in Figure [l^. A 
full calculation of the fields generated by the field coil 
would require a three-dimensional solution of coupled 
London's and Maxwell's equations in this geometry. To 
estimate the errors associated with approximating the 
actual field coil geometry by an infinitely thin circular 
wire, we consider the idealized geometry of Fig. [TJd: the 
field coil is assumed to be an incomplete, infinitely nar- 
row circle of radius a, which connects with infinitely long, 
infinitely narrow leads with spacing 5, and which carries 
a current /. For these calculations, we take a=8.4/im 
and s = 7.3/im. s was taken as the geometric mean of 
the outside 13/im and inside 1.2 jam widths of the leads 
in the susceptometer layout. 
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Using Biot-Savart: 



B 



jhqI C dl X f 
47r J \r\^ 



the contribution Bz^c^^ the z-component of the field from 
the incomplete loop is given by 



(A.l) 



fiol f'^^ a'^ — ay sin — ax COS 

'''^^Jo, {{x-a cos 0)^ -^{y-a sin 0)^ + ^2)3/2 ' 

where Oi — cos~^(5/2a). The contribution B^^i from the leads is 



/ip/ 2(g - 2y)(2{x - xq) + - 2yY + 4((x - xo)^ + z^)) 
27r ((5 - 2^)2 + 4^2) ^^(5 - 2^)2 + 4((x - Xo)^ + z^) 



z=0 



■ Circular loop 
Incomplete-circular loop 

- Leads 

- Incomplete-circular loop+leads 




FIG. 6: Calculated ^-component of the field Bz^ divided by 
the current / for the model of Fig. [TJ), for a circular loop 
of the same radius (dot-dashed curve), the contribution from 
the incomplete circular loop (dotted curve), the contribution 
from the leads (dashed curve) and the sum of the previous 
two (solid line). The x-axis is oriented along the leads, with 
the leads coming in along the positive x-direction. This cal- 
culation assumes the pickup loop and the field coil are in the 
same plane. The field at the center is 7.5% larger for the 
circular loop model than for the incomplete circle plus leads 
model. 



where = a cos d\ . Figure [6] plots the z-components 
of the fields from a circular loop (^i = 0), from an in- 
complete circle, from the leads, and the sum of the in- 
complete circle and leads, using the parameters above, 
and assuming z=0: the field coil and the pickup loop in 
the same plane. The field from the circular coil model 
is 7.5% higher than that from the incomplete circle plus 
leads model at z = 0. At a height of z = 3/im, more 
appropriate for calculating a susceptibility, the error is 
6.8%. 

To calculate the SQUID susceptibility using the 




FIG. 7: Calculated magnetic scalar potential (ps{r, 0), divided 
by the current through the loop, for the model of Fig. [l]: with 
a=8.4/im, s = 7Aiim and zq = 1.5/im. 



present formalism requires calculating the magnetic 
scalar source potential (/^^(r, 0) for the above geometry. 
Converting Eq. 36 of Reference! 13 1 to SI units: 



^s(^O) 



dp? 



{{r-r'Y^zlf/^ 



(A.4) 



This is difficult to integrate over an arbitrary geometry 
analytically. Instead we did the integrations and Fourier 
transforms numerically. Figure [7] shows the results for 
(Ps{f^O)/I (which is dimensionless) for a = 8.4/im, s = 
7. 3/im and zq = Lb/j^m. 

The source field is given by Hg = V(^s(r, z). Figure 
compares the results obtained using Biot-Savart with the 
gradient of the scalar potential of Fig. [T] 

The Fourier transform of the response field at z = 
is given by 



-kzo 



(A.5) 
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FIG. 8: z-component of the source field, divided by the cur- 
rent, with a=8.4/im, s — 7 Aiim and zq — 1.5/im, calculated 
using Biot-Savart vs taking the gradient of the scalar poten- 
tial. Also shown for comparison is a cross-section through the 
response field image of Fig. [5] 
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FIG. 9: ^-component of the response field, divided by the 
current, with a=8.4/im, s — 7Aiim and zo — 1.5/im, t=10/im, 
and A — O.l/xm. 



where (pr,i{k) is given by Eq. |5j Figure [9] shows the 
results for the response field using a=8.4/im, s = lAjim 
and zq = 1.5/im, a sample thickness oi t = 10/im, A = 
0.1/im, fi2//^o = 1, Ms/mo = 1* the strong diamagnetic 
shielding, bulk, non-paramagnetic limit of Eq. [5j Table 
\n\ shows some selected comparisons of calculations using 
the circular wire model with the incomplete circle plus 
leads model for the field coil. 

The first 5 columns of Table [TTl are the field coil radius 
a, the height of the susceptometer above the sample sur- 
face zo, the spacing between the leads s, the thickness 
of the paramagnetic superconductor t , and the London 
penetration depth A. Hzifi)/! is the field at the center of 
the coil divided by the current through the coil. This is 



approximated in the circular wire model for the field coil 
by l/2a. Hzr{zo)/I is the response field at the center of 
the field coil (and pickup loop), divided by the current 
through the field coil. Finally Hzr{zo)/ Hz{0) is equiv- 
alent to (j){z)/(j)s^ the ratio of the sample susceptibility 
to the self-inductance in the limit where the pickup loop 
radius b « a. In all cases /i2 = Ms = y^o iii Table [TH 

In the thin diamagnetic limit (first 3 rows of the ta- 
ble) ^{z)/^s ^ -(a/A)(l - 2(z/a)/Vl + 4(z/a)2). In 
the bulk diamagnetic limit (j){z)/(j)s ^ + 4(z + 

A)V«^)^/^. These values are entered in the last column 
of the table. A comparison of the last 2 columns of the ta- 
ble shows that if one normalizes by the self-susceptibility, 
the analytic limits derived above for the circular field coil 
model agree with the full incomplete circle with leads 
model to within 10%, independent of whether the cur- 
rent is localized at the very inside of the field coil or at 
the outside of the field coil, and presumably for any cur- 
rent distribution in between. 

A more rigorous solution of the problem would solve 
London's equations for the current distribution in the 
field coil following, e.g. Brandt and Clem,^ then use 
those results to find the scalar potentials for a set of equi- 
current-density paths, and add them up with suitable 
weightings. However, the results of Table [TT| indicate that 
the results of such a complex calculation would not differ 
from the infinitely narrow, circular field coil model by 
more than 10 percent. 



b. Approximating flux in pickup loop by field at center of 
field coil times an effective area 

A simplification used in this paper is to approximate 
the fiux through the pickup loop by the field at the center 
of the field coil times an effective area. More traditional^ 
is to model the pickup loop area as composed of a circle 
of radius 6, co-planar and concentric with the field coil. 
An additional pickup area from the leads is approximated 
by a square of width and length offset from the center 
of the pickup loop circle by a length /\w. The square 
area contributes one third of the fiux passing through 
it to the total pickup loop fiux. For these calculations, 
we take a=8.4/im, s = 5/im, 6=1.8/im, w=A.biim. and 
AK;=1.5/im. 

Numerical integration of the field from the field 
coil gives an integrated fiux through the pickup loop 
(690c&o/A) that is 11% larger than the fiux 5^(0)A 
(612<l>o/A), where Bz{0) is the field at the center of the 
field coil, and A = 7r6^ + k;^/3 is the effective area of the 
pickup loop in the circular field coil model (see Appendix 
Taj, and 16% larger (693<l>o/A vs. 583<l>o/A) in the in- 
complete circle plus leads model, assuming z = 0. This 
would be appropriate for calculating the self-inductance 
of the susceptometer, and could help to explain why it 
is necessary to use a somewhat larger eflFective area for 
the pickup loop {22 iim? rather than lljim?) than the 
Ketchen moder^ gives. More appropriate for estimat- 
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8.33 X 10" 
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2 


4.42 X 10"^ 


0.836 
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12 


1.5 


13 


10 


0.1 


3.71 X 10- 


-2 


4.17 X 10" 


2 


3.27 X 10"^ 


0.882 


0.902 


6 


1.5 


1.2 


10 


0.1 


7.41 X 10" 


-2 


8.33 X 10" 


2 


5.50 X 10"^ 


0.742 


0.687 



TABLE II: Some results from the evaluation of Eq.'s [A]4l |A.5| and [5] for various parameters. The first 3 rows are in the thin 
film diamagnetic limit, and the last 3 are in the strong bulk diamagnetic limit. 



ing uncertainties in susceptibility measurements would 
be taking the pickup loop spaced 2z from the field coil. 
Taking 2z = 3/im gives 523^o/A vs 511<l>o/A for the 
circular loop model, and 504<l>o/A vs 489<l>o/A for the 
incomplete circular loop plus leads model: This approx- 
imation leads to about a 7% contribution to the total 
error in the mutual inductance. 



c. Analytical approximations 

Table [T] displays analytical approximations to the full 
expression Eq. [7|in various limits. These analytical ap- 
proximations can be evaluated much faster than the full 
expression, and in some limits the temperature depen- 
dence of the penetration depth A or the permeabilty /i2 
can be inferred from the data without curve fitting, aside 
from a multiplicative constant, but care must be taken. 



a) 



Thin paramagnet 



Figures 10 and 11 show contour plots of the fractional er- 
ror Inexact -(/>analytical|/0exact aSSOCiatcd with USiug Cach 

of the 5 analytical expressions in Table [l[ 

For the ^-doped STO data of Fig. [3^ the best fit value 
for A/a = 113, and t/a - 7x10-"^, so the error associated 
with using the analytical expression (Table [l^. Fig. [lOb) 
is approximately 1%. For the LaNiOs data of Fig. ^ 
t/a = 2.4 X 10~^, and the error associated with using Ta- 
ble [iJd, Fig. \iO\p is approximately 0.2%. The systematic 
errors associated with using the analytical expressions for 
the 2-DEL data of Fig. [5] are also negligible. 



o 
N 




0.05 0.1 0.15 0.2 0.25 0.3 
t/a 

Thin diamagnet 




0.05 0.1 0.15 0.2 0.25 0.3 
t/a 



d. Uncertainties in parameter values 

The largest systematic errors in determining material 
parameters such as the penetration depth A and the 
permeability /i of a permeable superconductor are un- 
certainties in the parameters such as the height of the 
SQUID susceptometer zq above the sample surface, and 
the change in sensor height with applied voltage dz/dV. 
Figure [12] shows estimates for the uncertainties in the pa- 
rameters A, zo, and dz/dV from fits to the (5-doped STO 
data of Fig. [3^. The gray-scale images in this figure dis- 
play the error square sum = Xln(^(^) ~ ^fit(^))^ for a 
three dimensional volume in parameter space, projected 
onto the three 2-dimensional planes A — 2:0, zq — dz/dV ^ 



FIG. 10: (a) Fractional error exact ^analytical | / 0exact aSSO- 
ciated with using the thin paramagnetic (Table [Id) analytical 
expression instead of the exact expression Eq. 7| assuming 
A ^ 00 and /is = Mo • (b) Fractional error for the thin diamag- 
netic (Table |l^) analytical expression, assuming — — Mo 
and zo/a = 0.2. 



and A — dz/dV by taking the minimum value of 'E? along 
each projection axis. The other two parameters, a ver- 
tical shift Scj) and a linear slope ^unear = o^z^ were op- 
timized for each pixel in the 3-dimensional parameter 
space. One way to estimate the uncertainty in the pa- 
rameters is to determine the region in parameter space 
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FIG. 11: (a) Fractional error Inexact - (/>anaiyticai|/0exact asso- 
ciated with using the bulk paramagnetic analytic expression 
(Table assuming fi3 — jiQ and A ^ oo. (b) Fractional 
error for the bulk strong diamagnetic expression (Table 
assuming t ^ oo and 112 — IJis — IJio- (c) Fractional error 
for the the bulk strong diamagnetic expression (Table |l)l), 
assuming zo/a — 0.2, and 112 — IJis — IJio ■ 



where 'E? is less than twice its minimum value. The 
global minimum value for dz/dV (2.9/im/V) is consis- 
tent with our knowledge of the physical properties of our 



z-bender at low temperatures. Fig. [12] shows that the 
best fit value for A depends sensitively on zq. We esti- 
mate from our knowledge of the tip-sample geometry that 
the sensor height 1/im < 2^0 < 2.5/im, which implies that 
700/im < A < 1100/im. As can be seen from Table [l^, 
the SQUID susceptibility in the thin diamagnetic limit 
is proportional to a/A, and therefore a systematic error 
in a will result in a proportional error in A. We consider 
it unlikely that our estimate of a is incorrect by more 
than ±20%, and therefore assign a further systematic er- 
ror of ±20% to uncertainties in the effective sizes of the 
field coil and pickup loop. The Pearl length can be re- 
lated to the density of superconducting carriers through 
Us = m*//ioe^A, where e is the elementary charge. Using 
m*=1.25mp,2^ results in = 3.8±2.5-l.l x 10^^1/cm^ 



Fig. 13 displays the error square sum for a three 
dimensional volume (x2^, ^0, and dz/dV) in parame- 
ter space, projected onto the three 2-dimensional planes 
X2t — ^0, ^0 — dz/dV ^ and X2^ — dz/dV for fits to the 
LaNiOs data of Fig. ^jp. If we assume that the suscepti- 
bility in LaNiOs arises from isolated paramagnetic spins, 
we can estimate the 2D substrate spin density Ng by us- 
ing X2t = l^oNs{gi^B)^J{J ^ l)/3kBT]^The systematic 
uncertainty in X2^ should again be proportional to our 
uncertainty in a. Assuming a ±20% uncertainty in a, 
g = 2 and J = 1/2 leads to Ns - 6.4 ± 5.1 - 2.3 x 10^"^ 
cm~^: The diamagnetic signal in (5-doped STO is 4000 
times bigger than the paramagnetic signal in LaNiOs, 
but the calculated superconducting carrier density is 70 
times smaller than the calculated spin density. 



2. Statistical errors 

The solid symbols and lines overlaid on the gray scale 
images in Fig. 



12 and 13 represent the best fit values 



and 95% confidence intervals for the parameters using a 
statistical bootstrap analysis. Briefly, in this analysis a 
random sampling of the data was generated, with sub- 
stitutions, to produce the same number of points as the 
original set. This set was fit to the model allowing all 
5 parameters to vary, best fit parameters were recorded, 
and the procedure was repeated 5000 times. A histogram 
of the best fit parameters was generated, and confidence 
interval limits were set at the 2.5% and 97.5% levels. 
In the case of the J-doped STO data of Fig.'s [3k 



and 12 it appears that the statistical uncertainties are 
smaller than the uncertainties associated with our im- 
precise knowledge of the sensor height zq . For the case of 
LaNiOs of Fig.'s|3]3 and[T3]the bootstrap statistical anal- 
ysis indicates that the statistical uncertainties dominate, 
as might be expected from the noise in the data. 
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FIG. 12: Plot of the sum of the errors squared, divided by the 
global minimum (S^/S^i^)? for the ^-doped SrTiOs sample 
data of Figure [sj), fit to the thin diamagnetic expression (Ta- 
ble [Je)), projecting the minimum value in the A, zo, dz/dV 
parameter space, taken along the third axis, onto the A — zq 
plane (a), the A — dz/dV plane (b), and the zq — dz/dV 
plane (c), taking fixed values a=8.4/im and 6=2.7/im. The 
global best fit values are A = 954/im, zq = l.T/xm, and 
dz/dV = 2.8/im/y. The solid symbols and lines are the 
best fit and 95% confidence limits for the parameters from 
a statistical bootstrap analysis. 



FIG. 13: Plot of the sum of squares error, divided by the 
global minimum (S^/S^i^)? for fits of the LaNiOs data of 
Figure [sji to the thin paramagnetic limit expression (Table 
IJb)), projecting the minimum value in the X2t, zo, dz/dV 
parameter space, taken along the third axis, onto the X2t — zq 
plane (a), the x^t — dz/dV plane (b), and the zq — dz/dV 
plane (c), taking fixed values a=8.4/im and b=2.7/im. The 
global best fit values are x^t = 1.3 x 10~^/im, zo = l.Q/xm, 
and dz/dV = 2.6/im/V. The solid symbols and lines are most 
probable values and 95% confidence limits for the parameters 
from a statistical bootstrap analysis. 
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